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Abstract 

We give an algorithm write down all conformally invariant differential operators acting 
between scalar functions on Minkowski space. All operators of order k are nonlinear and are 
functions on a finite family of functionally independent invariant operators of order up to k. 
The independent differential operators of second order are three and we give an explicit real- 
ization of them. The applied technique is based on the jet bundle formalism, algebraization 
of the the differential operators, group action and dimensional reduction. As an illustration 
of this method we consider the simpler case of differential operators between analytic func- 
tions invariant under the modular group. We give a power series generating explicitly all the 
functionally independent invariant operators of an arbitrary order. 

1 Introduction 

It is well-known that the Maxwell equations are conformally invariant. This motivate an perma- 
nent interest in studying the conformal classes of metrics, conformally invariant operators and 
structures. The Maxwell differential operator is linear and so it is a splitting operator between 
two (infinitely dimensional) linear representation of the conformal group. There are many papers 
in the literature treating the splitting operators between different representations of the conformal 
group. These operators generalize the Maxwell operator with respect of this property of invariance. 
These considerations are based on the studying of the description and the structure of the linear 
representations and their subrepresentations. 

The aproach applied in the present work is different. It is based on jet bundles technique ( 3 ), 
(0) and (0). The jet lifting of functions (or sections) plays a role of an universal differential op- 
erator. The differential operators are viewed as a composition of a jet lifting and a fibre preserving 
map (a morphism for the case of linear operators) between vector bundles. This algebraise the 
differential operators and after this the invariance means invariance of these maps (morphisms). 
During the next step the technique of group action and dimensional reduction are used to describe 
the invariant fibre-preserving maps (sections) . A crucial point is the study of the action of station- 
ary subgroup of some point on the jets of smooth sections at this point. This contains elements of 
the Catastrophe theory in the sense of R. Tom and Arnold (see [T]and [2]). If the steps have been 
done explicitly this prescription gives all invariant operators including the nonlinear ones. 

In our case (differential operators acting on scalar functions on Minkowski space) the main 
result is the following: all invariant operators are nonlinear of order equal or bigger than two 
(we exclude the trivial case of zero order operators). For second order (n — 2) there are three 
independent differential operators fsee 1351 l3*fjl and l3T|) . Any other invariant operator of order two 
is a function of three variables of them. ASimilarly, for order n — 3 there are 23 independent 
operators — the previous three and 20 operators more of order three. Every invariant operator up 
to order three is a function of these universal operators. For the case of an arbitrary finite order 
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n there exist a finite complete system of independent conformally invariant operators up to order 
n i.e. they generate all other invariant operators up to the corresponding order. 

In the simpler case of differential operators between analytic functions, invariant under the 
modular group the result is similar. All invariant operators are nonlinear. There are n — 2 inde- 
pendent operators of the order up to n. All others of are functions of them. We give a power series 
generating the family of independent invariant differential operators. The invariant operator of 
order three is closely related to the Schwarz derivative (there is a difference because the Schwarz 
derivative takes values in the quadratic differentials). 

2 General Scheme 

Let £ = (E,p.M) be a vector bundle where E is a m + n-dimensional manifold and M is a Tri- 
dimensional manifold. Smooth sections M — > E build an infinite dimensional space denoted by 
G°°(£). We will often use adapted coordinates (x^,u a ), i.e. coordinates satisfying the relation 
p(x^, u a ) — (x^), [i = 1, . . . ,m, a = 1, . . . , n. 

A connected Lie group G acts on £ by bundle morphisms 

V g ■ C : I, , Diff(E) : V, :i, • &, (x ) , 

where £ x = p^ 1 (x) is the fibre over x G M and t g G Dif f(M) is the projection of a morphism T g . 
In adapted coordinates (x^, u a ) the action of G reads 

(^uVKWW (s)« 6 (s))- 

The group G has a natural action on G°°(£) 

G x G°°(£) 3 (g, V) — <?(V) e C°°(0. 

9{i , ){ x ) '■— Tg (tp (t~ 1 (a;))) . (1) 
In the local coordinates (a; M , u a ) this action is given by 

g(mx) = (T g T b (t- 1 (x))ip b (t^ix)). 

A natural problem is the description of the vector subspace C°°(£)g C C°°(£) of all G- invariant 
sections. The invariance condition g(ij)) — ip locally looks as 

{T g ) a b {t-\x))i> b {t-\x))=r{xl 

or cquivalently 

(T g ) a b (x) 1 p b (x)=r(t g (x)). (2) 

In general a description of all G-invariant sections is hardly possible but under some natural 
requirements imposed on the group action it may be achieved. There is a "smaller" reduced 
bundle £g that smooth sections (without any restriction) in it are one-to-one correspondent with 
G-invariant sections in £ (elements of G°°(£)g). The abstract algebraic construction of £g consists 
of two steps. Consider the stationary subgroup H XQ = {h € G\th(xo) = xq} of a point xq € M. 
For x = Xq the G-invariance condition © is 

(T h )l(x )i> h (x )=ip a (x ) , heH X0 . (3) 

As a matter of fact, this is a restriction on the values of G-invariant sections at any given point 
xq. Let s<(£);e = {1? G £,x \Th(u) — ~u,\/h G H Xo } be the vector subspace of all fixed vectors. 
We assume that the collection of all spaces st(£) x for all x G M is a vector bundle st(£) C £, called 
a stationary subbundle. This condition limits the action of G. Obviously, G°°(£)g G C°°(st(£)). 

Note: The explicit construction of the stationary subbundle st(£) is the crucial point where 
the new structure of the reduced bundle arises. This is the most difficult step in our approach. 
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The second step consists in taking the quotient of the base M for the bundle st(£). We suppose 
that the projection t g has uniform orbits in the base M and M itself represents a total space of a 
smooth locally trivial bundle (M, tt, M/G) of homogeneous spaces. It is another requirement on the 
action of the group G. Let's consider an orbit of G linking up x, y G M, i. e. 3 g G G : y = t g (x). If 
tp G C , °°(^)g the value ip{y) is uniquely determined by the value of the section at x (in accordance 
to the invariance condition 

r(y) = r(t g (x)) = (T g r b (x)^ b (x). (4) 

Let (M, ix, M/G) be a trivial bundle and N be a global section i.e. N C M is transversal to 
the orbits. Because of the relation Q, a G-invariant section ip is completely determined if we 
know the restriction iP\n(iP\n G C°°(si(£)|iv))- Moreover, if <p G C°°(si(£)|jv) then correctly 
induces an invariant section tp G C°°(£)g. Indeed, let y G M, then the orbit through y intersects 
the submanifold N only in one point x : y = t g (x) for some g G G. By definition, the induced 
G-invariant section is 

r( V ) := {T g ) a b (x) V \x). 

The element g G G is not uniquely determined by y and a; but since f(x) G st(£);c the value tp(y) 
doesn't depend on the specific choice of the group element. The restriction st(£)|jv is a coordinate 
realization of the bundle £g. If we consider another submanifold N' C M transversal to the orbits 
in M, the corresponding restriction st(£)|jv' is another coordinate realization. There is a canonical 
isomorphism st(£)jy « st(£,))N induced by the group action. This procedure sews the abstract 
reduced bundle from the coordinate realization. If (M, tt, M/G) isn't trivial the construction of £g is 
analogous but slightly complicated. We have to sew local coordinate realizations. Smooth sections 
in the reduced bundle are one-to-one correspondent with G-invariant sections. We shall transform 
the problem of describing invariant differential operators into a problem for characterization of 
invariant sections in appropriate jet bundles. 

Consider two bundles £ and r\ over the same base M. A group G acts on both of them by the 
same projection t in M. The action of G on G°°(£) and C°°(r]) induces an action on differential 
operators D : G°°(£) -> C°°{ri) 

g{D){i,):=g{D{g-\i ) ))). (5) 

We use for simplicity the same notation for the actions of G on £ and on rj . A differential operator 
is called G-invariant if it satisfies the following condition 

g(D)(iP)=D(g(i;)), V.g G G, V?/> G C°°(£). (6) 

The problem we consider is the description of all invariant differential operators. This problem 
can be reduced to the problem we studied before by using the jet bundle technique (for more 
details about jet bundles see |S1,E];E1> El)- Let a differential operator D : G°°(£) — > C°°(r]) be 
of order (up to) k and linear. We denote with J fc (£) the corresponding fe-jet bundle of £. For 
a local coordinate frame (x^,u a ) in £ there exists an induced coordinate frame in J fe (£) denoted 
by (a; M ,u a ,w^, ...,u^ lfl2 ...„ fc ), where the indices are ordered fi\ < fa < • • • < Mfe' Each vector 
u G J k (£,)x is a jet of particular section tp : M — > E, i.e. u = j k {ip) x so that 

u a = r(x), v* = d,r(x), . . . , < M2 ... Mfc = d^.^rix). 

Any linear differential operator is completely determined by its general symbol, i.e. the bundle 
morphism 3> : J k (£,) — > ?y (over the identity on M). Using some natural isomorphisms, one can view 
general symbols as sections in the tensor product (J fe (£)) <8> rj. The set of all linear differential 

operators of order up to k corresponds one-to-one to smooth sections C°° ((J k (Q)* ® vj ■ The jet 

lifting of the sections j k : C°°(£) — > C°° (J k (£,)) pl a Y s the role of an universal differential operator 
of order (up to) k. An arbitrary linear differential operator D : C°°(£) — > G°°(rj) of order k is a 
composition D = S> o j k . If D is a nonlinear operator then its general symbol is a fibre-preserving 
map 9 : J fc (£) -> r?. The action of the group G in the bundle £ induces another action of G in the 
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corresponding jet bundle J fc (£) (so-called jet lifting of the action). Thus we have an action of G in 
the tensor product (</ fe (£)) ® V- Then a differential operator is invariant if and only if its general 
symbol is an invariant section in (</ fe (£)) <g> r\- So the jet bundle technique gives an algebraization 
of the (linear) differential operators. If we construct the reduced bundle of (j fc (£))* ® rj this will 
provide a full description of all invariant operators. In the nonlinear case general symbol mathscrD 
is a fibre-preserving map. At any point x e M the restriction $) x : (J k (£,) x — » ?7x) is a smooth 
(nonlinear) map. The description of the G-invariant nonlinear operators is similar to the linear 

case. We study the action of the stationary subgroup H XQ on the space C°° (^(J k (0) Xo i^xo) — _ 

the space of all smooth maps J k {t;)xo — * Vx an d then we have to find the fixed elements in it. We 
use this scheme to describe all conformally invariant differential operators acting on Minkowski 
space. As a simple illustration we will consider the two dimensional algebraic conformal case. In 
both cases the base M is a homogeneous space. The reduced bundle £q consists of one fibre over 
one point. The crucial point is to find the stationary elements of C°°(^) Xo for only one point. 



3 Illustrative example 

We consider the space of analytic functions of a single complex variable and the differential op- 
erators between analytic functions. The analytic functions may be viewed as sections on trivial 
line bundles over the complex plane C, i.e. ( = (Cx C,p,C). The group GL(2, C) acts on C by 
rational transformations 

**(*)== zeC ' d) eGL ( 2 > C ) (7) 

GL(2, C) acts on the analytic functions by transforming the argument. In adapted coordinates 
(z, u) this action is 

9(z,u) = (^qr|> u ) ■ ( 8 ) 

The problem to solve is the description of all GL(2, C)-invariant differential operators. The invari- 
ance condition now reads 

D ('(=Ti))- w »(=ri)- v '"(" ^ CL(2X) ' 

The complex plane is a homogeneous space, because translations in C are a subgroup of GL (2, C) 

9= ( n ? J t 9 { , z ) = z + b. 



x 1 , 

We choose zo = 0. Its stationary subgroup H is defined as the subset 

H = ^heGL(2,C) h=(^ a c J ^ a^oj. (9) 

The fibre J fc (£)o is the set of all A; -jets taken at z = of analytic functions. The fc-jet of an 
analytic function is its Taylor expansion up to order k. If w — j k (f)o € J k (£,)o then 



w = J2j\ uizl - (10) 



i=i 

where 



Ul = —> 1 = 1,2,...,*. 
We have assumed above for simplicity that uq = since it doesn't lead to any loss of generality. 
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The prolonged action of Ho on the fibre J k (£)o is given by the definition 



Hw) =j k (f 



cz+l 



j k u)°i k 



nz 



cz + l 



, h€H . 



(11) 



The bar over the right hand side of indicates that the composition is truncated, i.e. all 
monomials of higher order than k have been ignored. The jet of a composition of functions is the 
composition of their jets. On the other hand, the truncated Taylor expansion has the form 



cz + 



A— 1 



-J = az J2(- CZ Y 



(12) 



m=0 



The transformation h{u\, . . . , Uk) = (w±, . . . , Wk) is determined by the equation 

fc - / fc-i \ 1 k ^ 



7WIZ 



(13) 



1=1 \ m=0 / 1=1 

For example, the transformation of the jet of fourth order is 



w 2 

W4 



U\CL 

u 2 a 2 



2u\ac 



= 113a — 6u2<i c + 6u\ac 



(14) 



114a 



12u 3 a 3 c + 36u 2 a 2 c 2 - 24u 1 ac 3 . 



Since we consider only scalar functions we have to describe maps J fc (£)o — > C invariant under 
the action of Ho. Let J k (£)o/Ho be the quotient space and it : J fc (C)o - * J k (Oo/Ho be the 
canonical projection on it. A map S>q : J k (£)o — > C is said to be i/o-invariant if and only if there 
exists another map : J k (£)o/Ho — > C such that the relation @ = 2> o tt holds. One may rewrite 
that requirement in terms of commutative diagrams as follows 



J k (Oo 



J k tt)o/H 



\ 



/ 



(15) 



C 

In this sense the canonical projection 7r is a universal TJo-invariant map. The components of 7r (in 
any coordinates in the quotient space) are invariant. Any ii/o-invariant map is a function of the 
components of 7r. 

To describe the quotient space means we must find a canonical representative in every orbit 
of Ho. The jets with u\ = is an invariant subspace. Considering the general case with u\ 7^ 
in each orbit there is an unique representative with u\ — 1 and u 2 = 0. The projection on this 
canonical representative is given by the element h € Ho with a — l/ui, c — u 2 /2u\ more precisely 




The coefficients wi are coordinates in the quotient space 

u 3 _ 3 / u 2 \ 

W3 ~(uif Hhv 

W4 = j — O g- + g , 

(ui) (ui) (ui) 
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u 5 m M 2M4 , ori (u 2 ) 2 u 3 45 (u 2 f , . 

?«5 = g — 10 g + 30 y g (19) 

(ui) (ui) (ui) 2 

The coefficients u>/ are the general symbols of the 7?o-invariant differential operators at z = 0. 
Since translations act in a trivial manner on the /c-jets they look exactly the same at any point of 
C. If we consider the infinite order jets by using the previous method we can obtain the following 
generating power series 

= z + Y i j^kZ k . (20) 

fc=3 

The functions wi — Wi(ui, u 2 , . . . , ui), I = 3, 4, . . . are components of the canonical projection it. 
According to the previous remark the corresponding differential operators are universal (a complete 
system of invariants). Any invariant differential operator up to order k is a function of them. The 
invariant differential operators corresponding to the symbols ((171) and 1)18(1 look as follows 




and 



(/'(*))' (/'(*))° (/'(z)) D 

The expression in the parenthesis in 1(21(1 is the well-known Schwarz derivative. Thus we have con- 
firmed the classic result that the Schwarz derivative is an invariant differential operator with respect 
to the action of fraction-linear transformations. The Schwarz derivative takes values in quadratic 
differentials. The difference between 1(21(1 and the Schwarz derivative is the factor {f'(z))~ 2 be- 
cause Di(f) takes values in scalar functions. Further more, we have obtained that there are k — 2 
functionally independent GL(2, C)-invariant differential operators of order I < k. Any invariant 
differential operator of the same order is a function of them. 

Each invariant differential operator is nonlinear and it has a particular domain, i.e. it exists 
the requirement f'(z) ^ 0. 

Note: The fact that the fc-jet j k (f)o starts with enables us to exclude the trivial case when 
the differential operators are operators of zero order, i.e. functions of /. It is clear that such 
operators are for sure GL(2, C)-invariant. Since we exclude the invariant subspace of jets with 
Mi = we obtain in general differential operators which involve division by a power of f'(z). 



4 Minkowski case 

Let us consider the Minkowski space ^# = (R 4 ,?7) with the pseudoeuclidean metric tensor r/ = 
diag(— 1, 1, 1, 1). The conformal group C(l,3) consists of all diffeomorphisms ip : M 4 — > M 4 pre- 
serving the conformal class of ry, i.e. 

where w(x) is a smooth function( the Greek indices take the values 0, 1, 2, 3). The conformal group 
C(l,3) is induced by the following transformations 

1. translations: x^ — » x^ + a^, a S R 4 ; 

2. rotations: x 11 — > k^x 1 ' , A^ry^A^ = 7y Q( g; 

3. dilatation: x^ — » Aa; M , A > 0; 

4. special conformal transformations: x^ —* [x^ + b^x 2 ) / (l + 2x.b + x 2 b 2 ) , b € M 4 . 
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Note: The special conformal transformations have a correct global definition in the compactificd 
Minkowski space. 

The conformal group acts on the space of the smooth scalar functions R 4 — > R by transforming 
the argument. We are looking for differential operator between scalar functions invariant under the 
action of the conformal group. We are going to follow the scheme demonstrated in the precedent 
section. 

The Minkowski space is a homogeneous space. The stationary subgroup H C C(l, 3) for x = 
is induced by the transformations 

f . x» A£x v ; 

2. x» -* Xx^ ; 

3. x» -» (x» + b^x 2 ) I (f + 2x.b + x 2 b 2 ) . 

We must consider the action of H a on the jets J k (R 4 ) Q = |j fe (/)o | / e C°° (R) 4 , /(0) = o| and 

find the _ff -invariant functions J k (R 4 ) — > R, i.e. we have to describe the canonical projection 

tt: J fc (R 4 ) ^ J fc ( R4 ) / ff °- 

The first nontrivial case is fc = 2. The 2-jet of a smooth function / is the Taylor polynimial 

f(.f)o = u a x a + ^u aia2 x a 'x a \ (22) 

where (u a ,u aia2 ), a\ < a.2 represent coordinates in the fibre J k (R 4 ) - 

Let ip : R 4 — > R 4 be a diffcomorphism with a fixed point x — then its second order jet at 
x = is 

(j 2 ( V )o)"=^x" + ^^^. (23) 

assuming that det(^) ^ (A£ iU2 is arbitrary). The action of a diffcomorphism p on the space 
J 2 (R 4 ) is given by 



AV\ AV2 „, I Ay 7/ 



(24) 



A special case is the action of the stationary subgroup H$. We have for the prolonged action of 
the dilatation 



as well as for rotations 

( u » \ ^ ( K u " 
The infinite jet of special conformal transformations (p reads 

oo 

(j °°(<p(x))\ T = (x» + b»x 2 ) ]T(-l) fe (26.x + b 2 x 2 )' 

specifically for fc = 2 



(25) 



(26) 



k=0 



Thus the action is 



(j 2 {ip{x))\ y = x» -2x"b.x + b'"x 2 . (27) 

(28) 



/ilM2 / \ "/J1/J2 111 "112 ^"■/X2 w All T ^"-p" 'IH1H2 



where & M = ^fc". 
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We will describe the quotient space J 2 (M 4 ) Q / Ho by choosing a canonical representative from 
each orbit of Hq. Let (u a , u aia2 ) be coordinates in J 2 (R 4 ) - 

At this step we assume that u 2 := rj^u^u 1 ' < 0. By dilatation choosing A = \ j \J (—u 2 ) we 
obtain another representative of the same equivalence class 



u a /\/-u 2 

^a-lOL^I ( U ) 



(29) 



By a special hyperbolic rotation 



[4 S$- v\ (1 + «q) /It 2 J ' (30) 

where we denote by if the space component (y 1 , v 2 , v 3 ^\ of the 4- vector v, we reach to another 
representative 

(1,0,0,0) \/ e \ (31) 

In the case of 2-jets looking like (|31[l the action l|28|> reads 

e — > e 
wo/3 — > W /3 - 2a/3 
Wij ™ > Wij - 2a (5ij. 

This formula enables us to choose = wqp/2 and therefore we get another "more canonical" 
representative 

/ (1,0,0,0) \ 



w 



mil 



Each equivalence class of this type of jets has got such a representative. The subgroup 0(3) C 
0(1, 3) is still acting on the jets 113211 preserving their form. The action on symmetric 3x3 matrices 
is w — > BwB T , £? € 0(3). That is why we can choose the transformation (a different boost) turning 
the matrix w into a diagonal form 

(1,0,0,0) \ / (1,0,0,0) 

( \ ^ / o \ 

\ w mn J J \ \ diag(M, A 2 , A 3 ) / 

The unordered triple of eigenvalues (Ax, A2, A3) are coordinates in the quotient space J 2 (M 4 )o/-ffo: 
so they describe completely this quotient space. As coordinate frame we may choose the elementary 
symmetric polynomials 

o"x = Ax + A2 + A3, (T2 = AxA2 + Ax A3 + A2A3, 03 = AXA2A3. 

but we prefer working with the following frame 

Si (Ax, A2, A3) = AX+A2+A3 

5 2 (Ax,A 2 ,A 3 ) = (Ax) 2 + (A 2 ) 2 + (A 3 ) 2 (33) 
S3 (Ax, A 2 , A 3 ) = (Ax) 3 + (A 2 ) 3 + (A 3 ) 3 . 

because of its convenient form 

S fc (Ai,A 2 ,A 3 ) ^Tr(w k ), k = 1,2,3. (34) 
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The quantities S\,S2 and S3 are coordinates in the equivalence class of the jet that we started 
from. To obtain them as explicit functions of the initial jet we must take the composition 



Tr(w) 
Tr(w) 3 



The functions 2>]~ — Tr (w k ) = @>k {u a ,u aia2 ) are the general symbols of the invariant operators 
at x = 0. Since the translations act trivially on the jets the differential operators look in the same 
way at any point of the Minkowski space. The final result is 

+2 (35) 



D 2 (f) = (*?°V A -2f?^ a ( ^/ +2 

v 2 / d«f ^/ a a/J / 



#7 A / d a f d?fd x f d».f\ d aP f d\/j,f 



(v/) 2 (v/) 2 ; (v/) 2 (v/) 



(V/) 2 (V/) 2 (V/) 2 



(36) 
(37) 



where 



n (f] = r *ex» P « d aP f d x ,f d pa f V 2 / jrfdPffrf&f d*pf gy/ 
3UJ (V/) 2 (V/) 2 (V/) 2 + (V./) 2 (V/) 2 (V./) 2 (V/) 2 (V/) 2 ' 

aa dPfd\fd»fd?f au d?fd\fdrfd°f 
v (v/) 2 (v/) 2 77 (V/) 2 (V/) 2 + 
s a /^/ a A /^/ 

2 (V/) 2 (V/) 2 (V.f) 2 ' 

(v/) 2 = v a0 d a fdfsf, v 2 / = ^a Q/3 / 

So we have just obtained the differential operators having considered functions with time-like 
gradients. According to the general scheme these differential operators are functionally independent 
and universal, i.e. they generate all conformally invariant operators of second order (defined on 
functions with time-like gradient). This procedure also contains an algorithm for calculating the 
higher order differential operators (defined on the same subset of functions). 

The natural framework of describing the conformally invariant differential operators involves 
the complexified and compactified Minkowski space. 

This technique is applicable to the case of n-dimcnsional pseudoeuclidean space too. There are 
no conformally invariant first order differential operators. The conformally invariant differential 
operators of second order are generated by n — 1 functionally independent differential invariants. 
For example, one of these invariant operators (involving functions with time-like gradient) is the 
following one 

L>(t) = 7T + (n — 2) 5 » a,p = l...n 

(V/) 2 (V./) 2 (V/) 2 
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